Fragmentation paths in dynamical models 
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We undertake a quantitative comparison of multi-fragmentation reactions, as modeled by two dif- 
ferent approaches: the Antisymmetrized Molecular Dynamics (AMD) and the momentum-dependent 
stochastic mean-field (SMF) model. Fragment observables and pre-equilibrium (nucleon and light 
cluster) emission are analyzed, in connection to the underlying compression-expansion dynamics in 
each model. Considering reactions between neutron-rich systems, observables related to the isotopic 
properties of emitted particles and fragments are also discussed, as a function of the parametrization 
employed for the isovector part of the nuclear interaction. We find that the reaction path, particu- 
larly the mechanism of fragmentation, is different in the two models and reflects on some properties 
of the reaction products, including their isospin content. This should be taken into account in the 
study of the density dependence of the symmetry energy from such collisions. 

PACS numbers; 25.70.Pq,24.10.Lx,24.10.Pa 



I. INTRODUCTION 

Due to the analogies between the nuclear forces and the 
Van-der-Waals interaction, the nuclear matter equation 
of state (EOS) foresees the possible occurrence of phase 
transitions from the liquid to the vapour phases [l|, 
Already in the 80's, these phenomena were linked to the 
experimentally observed multifragmentation mechanism, 
i.e. the sudden break up of an excited nuclear system 
into many pieces [1, 01 • However, in the search for sig- 
natures of a phase transition, one has to consider that 
nuclei are finite systems. New features may appear in 
addition to standard thermodynamics. Moreover, the nu- 
clear phase diagram is explored with the help of nuclear 
reactions, where out-of-equilibrium effects may also be 
present. These considerations have lead, from one side, 
to significant developments of the thermodynamics of fi- 
nite systems and to new definitions of the phase transi- 
tion, following the concepts of statistical physics [1, Q . 
On the other hand, many efforts have been devoted to the 
study of the reaction dynamics and to the characteriza- 
tion of the fragmentation mechanism, also in connection 
with the possible appearance of the signals expected for a 
phase transition at thermodynamical equilibrium (7l-[l3| . 

It is generally believed that, in a heavy ion collision 
at Fermi energies, due to the initial compression and/or 
thermal excitation, the composite nuclear system may ex- 
pand and reach density and temperature values inside the 
co-existence region of the nuclear matter phase diagram 
0, S [3 • In violent collisions multifragment configura- 
tions are characterized by a large degree of chaoticity and 
a huge amount of the available phase space is populated. 
Finally several features observed in the exit channel can 
be related to the thermodynamical (equilibrium) proper- 
ties of phase transitions in finite systems, independently 
of the specific mechanism that has driven the whole pro- 
cess However, there remain some aspects that 
refiect the nature of the reaction dynamics and of the 
fragmentation path, mostly related to out-of-equilibrium 



effects. Indeed the interplay between the compression- 
expansion dynamics and the onset (and nature) of clus- 
terization affects significantly the fragment kinematical 
properties and the appearance of collective flow effects. 
A thorough study of this dynamics also allows one to ac- 
cess information on the nuclear matter compressibility. 
Moreover, the specific features of the fragmenting source, 
such as its mass, charge and excitation energy, are also 
strongly depending on the pre-equilibrium dynamics and 
on the time instant where clusterization sets in. 

Hence, a quantitative understanding of multifragmen- 
tation data requires a careful investigation of the whole 
dynamical path. Due to the complexity of the nuclear 
many-body problem, two main lines of approximations 
have been followed so far. On one side, the class of molec- 
ular dynamics (MD) approaches employ, to represent a 
many-body state, a product of single-particle states, with 
or without antisymmetrization, where only the mean po- 
sitions and momenta are time dependent 'l5'-'20'|. In al- 
most all these approaches, the width is fixed and is the 
same for all wave packets. The use of localized wave 
packets induces many-body correlations (analogous to 
those in classical dynamics) in the particle propagation 
in the nuclear field, as well as in hard two body scat- 
tering, that is treated stochastically. On the other side, 
mean-field approaches (and stochastic extensions) follow 
the time evolution of the one-body distribution function 
(the semi-classical analog of the Wigner transform of 
the one-body density matrix), according to approximate 
equations where higher order correlations are neglected 
(mean-field approximation), apart from the correlations 
introduced by the residual two-body collisions. In the 
stochastic extension of the model, a fluctuation source 
term is added to the average collision integral, to ac- 
count for the stochastic nature of two-body scattering 

A comparison of the predictions given by models of 
the two classes, concerning multifragmentation scenar- 
ios at Fermi energies (30-50 MeV/nucleon), has been re- 



2 



cently undertaken [23|. The models considered are: the 
stochastic mean-field (SMF) model including momentum 
dependence [l^ and the antisymmetrized molecular dy- 
namics (AMD) model It was observed that, while in 
the SMF case fragment emission is more likely connected 
to the spinodal decomposition mechanism, i.e. to mean- 
field instabilities, in the AMD approach many-body cor- 
relations have a stronger impact on the fragmentation 
dynamics, leading to the earlier development of density 
and momentum fiuctuations. 

In the present paper we will discuss quantitatively the 
impact of the different approximations employed in the 
two models on several particle and fragment observables 
of experimental interest. This also allows one to establish 
up to which extent some of the observed features may be 
considered as more robust or general, i.e. not much de- 
pending on the details of the fragmentation path and of 
the models under consideration. By studying reactions 
with neutron-rich systems, we also investigate isospin ob- 
servables, in connection to the density dependence of the 
symmetry energy and the underlying reaction dynamics. 

II. DESCRIPTION OF THE MODELS 

A. Basic framework 

In both SMF and AMD approaches, the time evolution 
of the system is described in terms of the one-body dis- 
tribution function / (or a Slater determinant), as ruled 
by the nuclear mean-field (plus Coulomb interaction for 
protons) and the residual interaction, i.e. hard two-body 
scattering. The equation can be represented in the form 
of so-called Boltzmann-Langevin equation [21I, 

^^{H[f]J} + I[f] + 5I[f], (1) 

where the coordinates of spin, isospin and phase space 
are not shown for brevity. H[f] is the self-consistent one- 
body Hamiltonian, /[/] is the average two-body collision 
integral, and 5I[f] stands for the stochastic source term 
[21, 22]. 

The effect of the stochastic term 5I[f] is considered for 
the value of the distribution function in each phase 
space cell a. For example, one may decompose the phase 
space into square cells of the volume {2TTh)^. In the 
case of local equilibrium, the source term SI[f] induces 
the fluctuation of /„, the variance of which is given by 
(A/^) = fa{l ~ fa)- The correlations between different 
cells are usually assumed to be small, but minimal corre- 
lations are introduced for the conservation laws. In the 
approximate treatment of SMF presented in [2^ [23] , the 
fluctuations are projected and implemented only onto the 
ordinary space, agitating the spatial density proflle from 
time to time during the reaction after the local thermal 
equilibrium is reached. 

It should be noted that the variance (A/^) — fa{l — 
fa) with the average value fa is achieved by choosing 



fa + ^fa = 1 and with the probabilities fa and 1 — fa, 
respectively. Namely the fluctuation will choose A fully 
occupied phase-space cells, with A being the number of 
nucleons in the system, while the remaining part of phase 
space becomes empty. AMD represents this situation by 
using Gaussian wave packets oc e~^'^^'"~^') -(p-pO /2'i i' 
{i = 1, . . . , A). Each wave packet has the minimum un- 
certainty AxAp = ^fi and therefore can be regarded as 
a phase space cell. The width parameter v = (2.5 fm)~^ 
defines the shape of the cell in phase space. 

In AMD, the right hand side of Eq. ([T]) is usually de- 
composed in a different way for the convenience in writ- 
ing a stochastic equation for the wave packet centroids. 
The mean field effect is decomposed to the motion of the 
centroids and the change of the shape of the one-body 
distributions: {H,f} = {H, /}cent + {H, /} shape- The 
first term {H, /}cent is calculated by employing the fully 
antisymmetrized wave function [l7| . while {-ff, /}shapc is 
app roximately calculated by using a test particle method 
[281 [29} . The stochastic source term is also separated as 
SI = SIcoW + <^ /split ■ The part <5/coii is the stochastic ef- 
fect included in the stochastic collisions {I + SIcoii), that 
move the centroids of the wave packets by choosing the 
scattering angle randomly. When a nucleon collides with 
another, Slgpht randomly selects a Gaussian wave packet 
from the single-particle state, that has been changing its 
shape according to {H, /}shapc, by splitting it into possi- 
bilities of Gaussian wave packets (quantum branching or 
decoherence) 

Thus the main difference between the two models 
lies in the implementation of the stochastic term, i.e., 
whether the fluctuation in the momentum space is in- 
tegrated out (in SMF) or not (in AMD). Moreover, in 
AMD, when a collision happens, two entire nucleons are 
moved to new wave packets in phase space, while it is 
not necessarily the case in SMF. Therefore the fluctua- 
tion is expected to have stronger impact on the collision 
dynamics in AMD than in SMF. 



B. Effective interaction 

The isospin and momentum dependent effective inter- 
action employed in the SMF model is derived via an 
asymmetric extension of the GBD force [lO, [HI, lead- 
ing to the BGBD potential [H [llj. In the AMD 
model, we employ the Gogny interaction which are com- 
posed of finite-range two-body terms and a zero-range 
density-dependent term. 

The corresponding nuclear matter EOS can be writ- 
ten as {E/A){p,/3) = iE/A)ip) + Esyrr.ip)P^ + 0(13^) + .., 
where the variable /3 = (p„ — Pp)/p defines the isospin 
content, being pr (r — n, p) the neutron or proton density 
and p the total density. We use a soft equation of state 
for symmetric nuclear matter (compressibility modulus 
Knm{po) = 215 MeV for BGBD and K^m = 228 MeV 
for the Gogny force). In Fig. [T] we report the EOS, as 
a function of the density p, for symmetric nuclear mat- 
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FIG. 1: Equation of state (EOS) of symmetric nuclear mat- 
ter, corresponding to the Gogny (full) and BGBD (dashed) 
interactions. 



ter, /? = 0, at zero temperature. We can easily adjust 
the parameters of the interaction in order to change the 
density dependence of symmetry energy (Esym) without 
changing the EOS of symmetric nuclear matter. The 
density behaviour of the symmetry energy, Esymip), is 
shown in Fig. [2] We remind that this quantity gets a ki- 
netic contribution directly from basic Pauli correlations 
and a potential part from the highly controversial isospin 
dependence of the effective interactions [IJ, |35( . For the 
BGBD force, we adopt two different parametrizations of 
the symmetry energy: "asy-soft", that gives a flat be- 
haviour around normal density, followed by a decreasing 
trend at large density, and "asy-stiff" , where the symme- 
try energy exhibits an almost linear increase with den- 
sity The Gogny (Dl) force [s^l has a quite simi- 
lar Esymip) to the "asy-soft" parametrization of BGBD, 
while another parametrization [SS*] of the Gogny force 
has Esymip) similar to the "asy-stiff" version of BGBD. 
Hence, in the following we will use the labels "asy-soft" 
and "asy-stiff" for the Gogny forces as well. 

Within the considered form of the nuclear interaction, 
the isovector terms depend explicitely on the nucleon mo- 
mentum k = p/h, leading to the splitting of neutron 
and proton effective masses. The mean-field potential 
felt by neutrons and protons in asymmetric nuclear mat- 
ter (/? — 0.2) is presented in Fig. 3, as a function of 
the momentum k, for three density values and for the 
adopted interactions (asy-soft and asy-stiff). As an effect 
of the momentum dependence, the difference between 
neutron and proton potentials decreases with fc, becom- 
ing negative at high momenta, especially in the asy-soft 
case. This corresponds to the proton effective mass be- 
ing smaller than the neutron one in neutron-rich mat- 
ter. Moreover, one can see that the difference between 
neutron and proton potentials is larger, at low density, 
in the asy-soft case, corresponding to the larger value 
of the symmetry energy (see Fig. [5]) , while the opposite 




FIG. 2: Parametrizations of the density behaviour of the 
symmetry energy adopted in the calculations: Asy-soft (black 
lines), Asy-stiff (grey lines). Full lines correspond to the 
Gogny forces, while dashed lines are for the BGBD forces. 



holds above normal density. The results of the BGBD 
and Gogny forces look very close to each other below 
and around the Fermi momentum, for both neutrons and 
protons. However, due to the different shape of the mo- 
mentum dependence in the BGBD interaction, with re- 
spect to the Gogny force, at high momenta (fc > 3 fm^^) 
potential are less attractive in the BGBD case. This dis- 
crepancy should not affect our comparison since we are 
interested in reactions at Fermi energies. 



C. Two-nucleon collision cross sections 

As the two-nucleon collision cross sections icFpp = (Jnn 
and (Tpn), we use the energy- and angle-dependent values 
in the free space with the maximum cut-off of 150 mb 
in both SMF and AMD calculations. We have confirmed 
that the degree of stopping reached in the reaction is 
similar for both calculations [13]. However, it should be 
noticed that these cross sections are larger than those 
adopted by the AMD calculation in Ref. |2a |. 



D. Fragment identification 

In the SMF model, the reaction products are recon- 
structed by applying a coalescence procedure to the one- 
body density pir), i.e. connecting neighboring cells with 
density p > ^po ("liquid" phase). In this way one can 
also identify a "gas" phase (p < |po), associated with 
particles that leave rapidly the system (pre-equilibrium 
emission) and/or are evaporated. Once fragments are 
identified, from the knowledge of the one-body distri- 
bution function, it is possible to calculate their mass. 
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FIG. 3: Mean-field potential, for neutrons (grey lines) and 
protons (black lines) in asymmetric matter (/3 = 0.2) at den- 
sity p = 0.085 fm"^ (panels (a) and (b)), p = 0.17 fm"^ 
(panels (c) and (d)), p — 0.34 fm~^ (panels (e) and (f)), as 
a function of the momentum k. Full lines correspond to the 
Gogny forces while dashed lines are for the BGBD forces. Left 
panels: asy-soft interaction. Right panels: asy-stifT. 



charge and kinematical properties. This provides rea- 
sonable results for the description of the ground state of 
medium-heavy nuclei and for excited primary fragment 
properties. Fragment excitation energies are calculated 
by subtracting the Fermi motion, evaluated in the local 
density approximation, from the fragment kinetic energy 
(taken in the fragment reference frame) 0, [IB] • 

In the AMD model, for the state at a given time of 
the reaction, fragments are recognized by connecting two 
wave packets if the spatial distance between their cen- 
troids is less than dc\. When fragments are recognized at 
a late time, such as i = 200 fm/c, we take dd — 5 fm, 
though the result does not depend on the choice except 
for extreme choices. When fragments are recognized at 
early times in order to characterize the dynamical evolu- 
tion of the reaction, different choices of dc\ define different 
observables. We take a relatively small value of del = 3 
fm in Figs. [HI El HH and [131 which allows an earlier recog- 
nition of particle emissions. 

Alternatively, we can employ the same fragment recog- 
nition method that is applied to the SMF calculation by 
applying a coalescence procedure to the density distri- 
bution p(r) calculated for the AMD wave function. We 
have confirmed that the result well agrees with the usual 
method with d^i = 5 f m for fragment charge distributions 
(i.e., for the separation of fragments), while the choice of 
del = 3 fm better agrees with the coalescence procedure 
for the light particle emissions (i.e., for the separation of 



liquid" and "gas") at early times. 



III. RESULTS 

We discuss some features of the fragmentation path 
followed in violent collisions at Fermi energies, as pre- 
dicted by the SMF and the AMD models. To investigate 
isospin effects and the sensitivity of the results to the 
parametrization adopted for the symmetry energy, we 
will simulate central collisions of a neutron-poor and a 
neutron-rich system: ^^^Sn -I- ^^^Sn and ^^^Sn -I- ^^''Sn, 
at 50 MeV/nucleon. Within our study of fragmentation 
reactions, we expect to test the low-density behaviour of 
the symmetry energy. As we will see in the following, 
the isospin degree of freedom can also be used as a good 
tracer of the reaction dynamics. An ensemble of about 
100 trajectories has been collected for both model calcu- 
lations. The impact parameter has been set equal to 0.5 
fm and the initial distance between the two nuclei (at the 
time i = 0) is 15 fm. Density contour plots in the reac- 
tion plane, as obtained in the two models for one event 
of the ^^^Sn 4- ^^^Sn reaction, are shown in Fig. [Hat sev- 
eral time steps. As one can see from the figure, for this 
kind of reactions, both models predict that the system is 
initially compressed. Then expansion follows and multi- 
fragment breakup is observed. One can also observe that 
the expansion is faster in AMD and that the amount of 
emitted nucleons (the "gas" phase) is larger in SMF. 



A. Compression-expansion dynamics 

In this Section we summarize the main findings of 
Ref. ^2j], concerning the trajectory followed by the sys- 
tem in the early stage of the dynamics, until fragmen- 
tation is observed. The reaction path can be character- 
ized with the help of one-body observables, such as the 
radial density profile and the radial collective momen- 
tum. The radial density at a given distance r is obtained 
by averaging the local density p(r) over the surface of a 
sphere of radius r. The radial collective momentum is the 
projection of the collective momentum at the position r 
along the radial direction, averaged over the surface of 
the sphere of radius r. These quantities are further av- 
eraged over the ensemble events. 

In the SMF calculations (see Fig. 5(a)) the behaviour 
of the radial density profile indicates that, after an ini- 
tial compression (t = 40 fm/c), the system expands and 
finally it gets rather dilute, due to the occurrence of a 
monopole expansion, generated by the initial compres- 
sion. The matter appears mostly concentrated within 
a given interval of the radial distance (see for instance 
the results at t = 100 fm/c), indicating the formation of 
bubble- like configurations (see also Fig. 2]), where frag- 
ments will appear, and the central region of the system 
is rapidly depleted. 
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FIG. 4: Contour plots of the density projected on the reaction plane calculated with SMF (top) and AMD (bottom) for the 
central reaction ^^^Sn + ^^^Sn at 50 MeV/nucleon, at several times (fm/c). The lines are drawn at projected densities beginning 
at 0.07 fm~^ and increasing by 0.1 fm~^. The size of each box is 40 fm. 
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FIG. 5: Density profiles, at several times, as obtained in the 
SMF (a) and AMD (b) case, for the reaction "2Sn+ "^Sn. 



In AMD calculations, the density profile in Fig. 5(b) 
shows the time evolution of the compression and expan- 
sion which is qualitatively similar to the SMF case. How- 
ever, we notice that AMD shows broader average density 
distribution than SMF as the system expands, pointing 
to a faster expansion in the AMD, as indicated qualita- 
tively by Fig. 21 This is confirmed also by the behavior 
of the collective momentum. An shown in Ref. [23|, in 
the SMF case, after the development of an almost self- 
similar radial flow, the collective momentum decreases 
again at the surface of the system, indicating the occur- 
rence of a counter-streaming flow, from the surface to- 
wards the interior, trying to recompact the system [39| . 
Correspondingly a bump is obtained in the radial den- 
sity distribution (see Fig. 5(a) at t = 100 fm/c). On the 
other hand, the collective momentum at the surface of the 
expanding system keeps almost unchanged in the AMD 
calculations [24] . The absence of deceleration effects may 
indicate that the system ceases to behave as homoge- 
neous while it expands, corresponding to the situation, 
of rather low average density, in which fragments have al- 
ready appeared and are distributed widely in space. This 
scenario is supported by the analysis of the density fluc- 
tuation variance, that reveals an earlier growth of density 
fluctuations (leading to an earlier appearance of clusters) 
in the AMD case [2J| . This can be interpreted as the rea- 
son of the faster expansion in AMD as observed in Fig. 0] 
at later times. We notice that early fragment formation 
is observed also in other N-body treatments, belonging 
to the class of molecular dynamics models, as shown by 
the quasi-classical calculations performed in Ref. [i^l . 

In conclusion, while in the SMF case the system spends 
more time, as a nearly homogeneous source, at low den- 
sity and fragments are formed most likely through the 
spinodal mechanism [2|, in the AMD case clustering ef- 
fects may be present at an earlier stage. This influences 
significantly also the amount of particles emitted prior 
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FIG. 6: Time evolution of the total number of neutrons (black 
lines) and protons (grey lines) contained in emitted particles 
with 1 < ^ < 4, in the case of the reaction ^'^^Sn + '^^''Sn. 
Dashed line: asy-soft parametrization. Full line: asy-stiff 
parametrization. The dot-dashed lines represent the time 
evolution of the total number of neutrons and protons con- 
tained in emitted particles with 1 < j4 < 15, for the asy-soft 
parametrization. Left panel: SMF results. Right panel: AMD 
results. 
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FIG. 7: The same as in Fig. [6] for the system "^Sn + "^Sn 



to fragment formation (pre-equilibrim emission), as dis- 
cussed in the following. 



B. Time evolution of nucleon and cluster emission 

During the first stage of the reaction, hard two-body 
scattering plays an essential role and pre-equilibrium 
emission is observed, i.e. nucleons and light particles are 
promptly emitted from the system. This stage influences 
significantly the following evolution of the collision. In 
fact, the amount of particles and energy removed from 



the system affects the properties of the composite source 
that eventually breaks up into pieces. Hence, when dis- 
cussing multifragmentation mechanisms, a detailed anal- 
ysis of this early emission is in order. Figs. [5] and [7] show 
the time evolution of the total number of neutrons and 
protons contained in emitted nucleons and light parti- 
cles, with mass number < 4, as obtained in the two 
models, with the two asy-EOS considered, and for the 
two reactions. It should be noticed that this ensemble is 
mainly composed of unbound nucleons in the SMF case. 
As shown in these flgures, these particles leave the sys- 
tem mostly in the time interval between « 70 fm/c and 
« 120 fm/c. At later times the emission rate is reduced 
(see the change of slope in the lines of Figs. [6] and [71 As 
already pointed out in Ref. [131 j ^ striking difference be- 
tween the two models concerns the amount of particles 
emitted, that is larger in the SMF case. However, the 
average kinetic energy of this emission is similar, being 
20.72 MeV/nucleon in SMF and 21.95 MeV/nucleon in 
AMD. The difference observed in the two models could 
be connected to the fact that clustering effects and many- 
body correlations are more efficient in AMD, due to the 
nucleon localization, reducing the amount of mass that 
goes into free nucleons and very light clusters. This effect 
is connected also to the different compression-expansion 
dynamics in the two models, as discussed in the previous 
section. 

On the other hand, light IMF's, with mass number 
5 < A < 15, are more abundant in AMD. This is ob- 
served in Fig. [6l where the total amount of neutrons and 
protons contained in emitted particles with A < 15 is 
also displayed as a function of time, in the case of the 
asy-soft interaction (see the dot-dashed lines). In the 
AMD case, the emission of light IMF's starts already at 
around 70 fm/c and, at the time t = 250 fm/c, it repre- 
sents a noticeable fraction of the particles emitted in the 
considered mass range {A < 15), see the difference be- 
tween dot-dashed and dashed lines. Hence fragments are 
formed on shorter time scales in AMD, on about equal 
footing as light-particle pre-equilibrium emission, while 
in SMF this light IMF emission sets in at later times and 
is much reduced. However, one can see that the total 
amount of emitted nucleons belonging to particles with 
A < 15 (including free nucleons) is close in the two mod- 
els (compare the dot-dashed lines at the final time in the 
left and right panels of Fig. Hence one expects to see 
a similar production of IMF's with charge A > 15, as it 
will be discussed in the following. 

As general features, when comparing the two reactions 
(Fig. inland Fig. [7]), it is seen that neutron (proton) emis- 
sion is more abundant in the neutron-rich (poor) systems. 
Moreover, Figs. |S1 and [71 also show that a larger (smaller) 
number of neutrons (protons) is emitted in the asy-soft 
case, as compared with the asy-stiff case, corresponding 
to a larger repulsion of the symmetry potential for the 
soft parametrization. This can be taken as an indica- 
tion of the fact that pre-equilibrium particles are mostly 
emitted from regions that are below normal density (i.e. 
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FIG. 8: Charge distribution as obtained in AMD (full his- 
togram) and in SMF (dashed histogram), for the reaction 
^^"Sn+^^^'Sn, at t = 200 fm/c (left) and t = 300 fm/c (right). 
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FIG. 9: Ratio Rz between the fragment yield observed in the 
angular domain 60 < S < 120 (rescaled by a factor 2) and the 
total yield, as a function of the fragment charge, evaluated at 
t — 300 fm/c for the same system of Fig.[Hl Notations are the 
same as in Fig. [S] 



after t « 70 fm/c, during the expansion phase), where the 
symmetry energy is higher in the soft case (see Fig. 2). 



C. Fragment properties 

Now we move to discuss and compare the properties 
of primary fragments, as obtained in the two models. 
In Fig. |8] we present the charge distribution of primary 
IMF's {Z > 2) at two time instants: t = 200 and 300 
fm/c. Some of the fragments identified by the clus- 
tering procedure a.t t = 200 fm/c have exotic (elongated) 



shapes and they break-up into pieces at later times. This 
effect is present in both models and explains the dif- 
ference observed between the charge distributions cal- 
culated at t = 200 and t = 300 fm/c. In the SMF 
case the maximum yield corresponds to fragments having 
charge around Z — 10, as expected within a fragmenta- 
tion scenario associated with spinodal decomposition [3] . 
The fragmentation mechanism may be connected also to 
the small yield of the light IMF's that, as also discussed 
in Section III.B, is lower in SMF, compared to AMD. 
Hence the total mass that belongs to the "liquid" phase 
(IMF's) is reduced in SMF. This is the counterpart of 
the more abundant pre-equilibrium emission. The pro- 
duction of sizeable clusters {Z > 6) looks closer in the 
two calculations. This is true especially a.t t = 300 fm/c, 
where the tail at large Z, that is observed at t = 200 
fm/c and is more pronounced in AMD calculations, dis- 
appears. However, the emission of fragments with charge 
around 10 is slightly larger in SMF, while bigger frag- 
ments (with charge around 15-20) are more abundant in 
AMD. It should be noticed that these differences may 
be smoothened by the secondary decay process, see the 
results of Ref.[4l| for SMF and Ref.fH for AMD. 

Fragments are nearly isotropically distributed in space, 
in both models. This is shown in Fig. 9, where we dis- 
play the ratio, Rz, between the fragment yield observed 
in the angular domain 60 < 6* < 120 (rescaled by a fac- 
tor 2) and the total yield, as a function of the fragment 
charge. Apart from fluctuations of statistical origin, we 
notice a lack of small fragments in the selected angular 
domain in the SMF case, pointing to a slightly lower de- 
gree of stopping, compared to AMD, as also discussed in 
Ref. [131. We would like to stress that, especially in the 
case of SMF calculations, the degree of stopping reached 
in the collision is crucial in determining the following re- 
action path. Indeed, if the reaction time becomes too 
short, spinodal instabilities would not have enough time 
to develop Q. 

The larger total mass of the "liquid" phase obtained 
in AMD is reflected also in the IMF multiplicity, that is 
shown on Fig. [TOl In fact, the multiplicity of IMF's with 
charge Z > 2 is larger in AMD. However, if one selects 
only sizeable fragments, with charge Z > 6, closer mul- 
tiplicities are obtained in the two models (dashed lines), 
but still higher (by about 1 unit) in the AMD case. This 
difference may be expected from the fact that the total 
number of nucleons contained in these IMF's in AMD is 
still larger by about 15 than in SMF (compare dot-dashed 
lines in Fig. 6). 

As discussed above, the average kinetic energy of nucle- 
ons emitted at the pre-equilibrium stage of the reaction 
is nearly the same in the two codes. However, the larger 
amount of nucleons emitted in SMF causes a consider- 
able reduction of the energy available for the remaining 
fragmenting system. As a consequence, the energy (ki- 
netic -I- internal) stored into IMF's results to be larger 
in AMD. More precisely, the fragment intrinsic excita- 
tion energy is close in the two models and amount to 
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FIG. 10: IMF multiplicity distribution as obtained in the two 
models for the same reaction of Fig. [S] at t = 300 fm/c. Full 
lines: IMF's with charge Z > 2 are considered. Dashed lines: 
IMF's with charge Z > 6 are considered. 
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FIG. 11: IMF average kinetic energy as a function of the 
charge Z, as obtained in the two models at t = 200 fm/c 
(dashed lines) and at t = 600 fm/c (full lines), for the same 
reaction of Fig. [8] The solid and hatched grey areas mark the 
change during this time interval for AMD and SMF, respec- 
tively. 



e/^ w 3 ± 1 MeV (corresponding to temperatures T « 5 
MeV), while according also to the different compression- 
expansion dynamics (see Section III. A), the fragments 
have higher collective velocities in AMD. 

The fragment average kinetic energy is represented as 
a function of the fragment charge in Fig. [TT] The dashed 
lines correspond to the time instant i = 200 fm/c, while 
full lines are for fragments at t = 600 fm/c, where they 
have been accelerated by the Coulomb repulsion. One 
can notice that in the AMD case, these Coulomb acceler- 
ation effects are not so pronounced, since fragments are 



rather distant from each other already at < = 200 fm/c 
(see Fig. 2]). The shape of the average kinetic energy, with 
an almost linear increase up to charges around Z = 15, 
denotes the presence of a collective flow velocity, as also 
experimentally observed in similar reactions [42i] . Frag- 
ments with larger size, that probably emerge from the 
coalescence of two or more smaller objects and are more 
likely located closer to the system center of mass, are 
slowed down. The difference between the final average 
kinetic energy per nucleon in the two models amounts to 
20%. 



D. Isospin effects 

It is very interesting to investigate also the isotopic 
content of pre-equilibrium and fragment emission, in con- 
nection with the asy-EOS employed and the details of the 
dynamical models under investigation. Differences be- 
tween the two models are seen not only for the abundance 
of the pre-equilibrium emission, but also for its isotopic 
content. This is illustrated in Fig. [TH that shows the 
time evolution of the {N/Z)^a_s ratio of emitted particles 
with ^ < 4, in both models, for the two iso-EOS's and 
the two reactions considered. The different calculations 
present some common features. In fact, we observe in all 
cases that the two curves of different asy-EOS's cross at 
times t around 70-100 fm/c, which is more evident for 
the neutron-rich system. According to the compression- 
expansion dynamics followed by the composite nuclear 
system, the crossing is connected to the fact that the sys- 
tem gradually evolves from a compact shape (compres- 
sion, high density) where the symmetry energy is larger 
and the system likes to emit more neutrons in the asy- 
stiff case, to a dilute configuration (expansion, low den- 
sity) where the symmetry energy is larger in the asy-soft 
case. 

The results of Fig. [12] can be connected to the differ- 
ences observed along the compression-expansion dynam- 
ics in the two models. In fact, as discussed in Section 
III. A, within the SMF trajectory the composite nuclear 
system may enter lower density regions as a nearly ho- 
mogeneous source and for a longer time, compared to the 
path followed in AMD. Hence the moderate {N /Z)^a.s of 
the pre-equilibrium emission observed in SMF can be ex- 
plained by the rather small value of the symmetry energy 
at lower density. This explanation is also consistent with 
the more abundant emission obtained in SMF calcula- 
tions, since particles are less bound at low density. On 
the other hand, the (iV/Z)gas ratio in the AMD case is 
much larger, especially for the neutron-rich ^^'*Sn reac- 
tion. Due to the efficient clustering effects in AMD, pro- 
tons are trapped in fragments that appear at early times, 
optimizing the symmetry energy of the whole system. 

The same isospin effects can be also observed by look- 
ing at the isotopic content of fragments. Figure [T51 shows 
the (7V/Z)iiq content of the "liquid" phase (associated 
with the composite nuclear source at early times and with 
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FIG. 12: Time evolution of the N jZ content of the pre- 
equihbrium emission, as obtained in the two models (left and 
right panels), for the "^Sn+^^^Sn and ^^''Sn+^^**Sn reactions. 
Full lines: asy-stiff interaction. Dashed lines: asy-soft. 
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FIG. 13: Time evolution of the isotopic content of the "liquid" 
phase, as obtained in the two models (left and right panels), 
for the "^Sn + "^g^ ^nd ^^''Sn + ^^^Sn reactions. Full lines: 
asy-stiff interaction. Dashed lines: asy-soft. 



IMF's at later times), as obtained in the two models and 
with the two parametrizations of the symmetry energy. 
As expected already from the results concerning nucleon 
and light cluster emission, the "liquid" is more neutron- 
rich in SMF. One can notice that, for the neutron-poor 
system in Fig. [131 the (A^/Z)uq ratio starts to increase 
at early times, which is due to the repulsive effect of the 
Coulomb interaction and is related to the observation in 
Fig.[l2]that {NjZ)^^ is smaller than the NjZ ratio of the 
total system. This increasing trend of (iV/Z)iiq is more 
pronounced in the SMF model, due to the more abun- 
dant emission in this case. For the neutron-rich system. 



the initial increase of (iV/Z)iiq is weak and observed only 
in the SMF case. The (iV/Z)iiq ratio turns to decrease 
approximately when fragments appear and neutrons are 
emitted, lowering the fragment symmetry energy (isospin 
distillation |34i| ) . This effect is more pronounced in the 
soft case and happens at earlier times in AMD. Moreover, 
isospin distillation is stronger in the neutron-rich system. 

At the end, one can see from Figs. [1^ and [T^] that 
the difference between the predictions of the two models 
is larger than the difference between the results associ- 
ated with the two iso-parametrizations. This observation 
should be taken as a warning that the symmetry energy 
cannot be extracted from merely the NjZ content of the 
reaction products in a model independent way. 



IV. CONCLUSIONS 

In this paper we have undertaken a quantitative com- 
parison of the features observed in heavy-ion fragmen- 
tation reactions at Fermi energies, as predicted by two 
transport models: SMF and AMD. As far as observ- 
ables of experimental interest are concerned, one signifi- 
cant discrepancy between the two models is connected to 
the amount of pre-equilibrium emission, i.e. the energetic 
particles that leave the system at the early stage of the 
reaction. In AMD clustering effects appear to be more 
relevant, reducing the amount of free nucleons emitted, 
compared to SMF, in favour of a richer production of pri- 
mary light IMF's. The yield of sizeable primary IMF's 
(Z > 6) is close in the two models. However in SMF frag- 
ments with charge around 10 are slightly more abundant, 
while in AMD the tail at larger Z (around 20) is more 
pronounced. The shape of the SMF charge distribution is 
closer to the expectations of spinodal decomposition 
However, it should be noticed that these differences can 
be smoothened by secondary decay effects [4l[. In fact, 
both models are able to fit (final) experimental charge 
distributions (of IMF's) reasonably well [28,, i^]- In SMF 
fragment kinetic energies are smaller, compared to AMD, 
by about 20%. These observations corroborate the sce- 
nario of a faster fragmentation process in AMD, while in 
SMF the system spends a longer time as a nearly homo- 
geneous source at low density, emitting a larger amount 
of nucleons prior to fragment formation. 

Another interesting feature is that the N/Z ratio of the 
pre-equilibrium emission and of primary IMF's is differ- 
ent in the two models. For a fixed parametrization of 
the symmetry potential, the isotopic content of the emit- 
ted particles is systematically lower in SMF, where the 
emission is more abundant. As a consequence, IMF's are 
more neutron-rich in SMF than in AMD. The latter ob- 
servation leads to the conclusion that isotopic properties 
are largely affected by the reaction dynamics. This can 
be expected just from the fact that the symmetry energy 
is density dependent, so isospin observables should keep 
the fingerprints of the density regions spanned in the col- 
lision. However, the impact of the reaction path on these 
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observables may be more intricate. From this point of 
view, isospin properties can be also considered as a good 
tracer of the reaction mechanism and may contribute to 
probe the corresponding fragmentation path. Hence, the 
discrepancy between the two models may be ascribed es- 
sentially to the different compression-expansion dynam- 
ics and clusterization effects that, as shown by our re- 
sults, affect isoscalar as well as isovector properties of 
the reaction products. 

Therefore, the simultaneous analysis of several experi- 
mental observables in nuclear reactions at Fermi energies 
should help to shed light on the underlying fragmenta- 
tion mechanism and the corresponding role of mean-field 
and many-body correlations. In other words, a check of 
the global reaction dynamics could be a way to test the 
validity of the approximations employed in the dynami- 
cal models devised to deal with the complex many-body 
problem. One may look at suitable isoscalar observables 
such as fragment and particle yields and energy spectra, 
as well as at the isotopic content of the reaction prod- 
ucts. As shown in Fig. 6, the sharing between the yields 
of light IMF's and light particles {A < 4) is rather differ- 
ent in the two models, the latter being more abundant 
and less neutron-rich in SMF. In particular, the study 
of pre-equilibrium emission, which is not influenced by 
secondary decay effects, could really help to probe the 
reaction dynamics, since the difference between AMD 
and SMF predictions is rather large, see the results of 
Fig.s 6,7,12. Only when the model reliability has been 
established, one may undertake a deeper investigation 



of isospin observables, from which more detailed infor- 
mation on the symmetry energy and its density depen- 
dence can be accessed. Indeed, as shown in Fig.s 12-13, 
the differences between predictions corresponding to two 
commonly employed symmetry energy parametrizations 
are smaller than the differences associated with the two 
models considered here. 

Finally, it would be appealing to extend our com- 
parison of multi-fragmentation reactions to other ap- 
proaches introduced to follow the many-body dynamics 
fisl [Tsl . [20j . In such a context, new developments of 
stochastic mean field models, in the direction of intro- 
ducing fluctuations in full phase space [23i| and enhancing 
the role of correlations, are also of interest. 
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